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^v^j . The effect of thermal fluctuations in the resonance fluorescence of a three-level system is 

studied. The damped three-level system is driven by two strong incident classical fields near 
resonances frequencies. The simulation of a thermal bath is obtained with a large system 
of harmonic oscillators that represent the normal modes of the thermal radiation field. The 
time evolution of the fluorescent light intensities are obtained solving by a iterative method the 



Or 
O 



Heisenberg equations of motion in the integral form. The results show that the time development 
of the intensity of the fluorescence light is strongly affected by the interaction of the system 
with the thermal bath. 
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In the last fifteen years, quantum jumps in three-level systems has been extensively studied [Q-|13 . The three- 



> 

level system is a versatile model that has been used to study, among others, the correlations in emission of photons 
Q, resonant fluorescence |j, operation of a two-mode laser |Q, coherent pump dynamics squeezing properties 
I/"-) , of electromagnetic field p0[ , electron shelving ]l2|| , quantum measurements theory The interest on processes 

involving atoms with few energy levels recently increased even more with the possibility of to study experimentally 
| non-linear processes with ion trapping. 

In this paper we consider the fluorescence resonance of a three-level system with a coherent interaction with two 
driving fields and a incoherent interaction with a thermal reservoir. When the driving fields are turned on, the system 
is driven to a new non-equilibrium steady state. If the driven fields are turned off, the system returns to the original 
equilibrium state with the thermal bath. We assume that the atom is in a cavity where there is a quantized radiation 
field in thermal equilibrium with the atom and the cavity walls at a certain temperature. In our model of a three-level 
atom, the allowed transitions are only between the levels 1 and 2 and 1 and 3. The three-level system of this kind 
is known as the V configuration. We study the system driven by the interaction with two electromagnetic fields of 
frequencies u>u e u>2h near, respectively, the lo\ e u)% frequencies characteristics of the system. The three- level system 
also interacts with a radiation field, with which one maintains in thermal equilibrium. We suppose that the normal 
modes of the radiation field constitutes the thermal bath at a certain temperature. We obtain the equations of motion 
' of the dynamical operators in the Heisenberg formalism in the general form 

^ = -/3x + A(t) (1) 

that may be recognized as equations that describes a damped system subject to fluctuating forces, accordingly 
Langevin's theory of Brownian motion p2],p3|. In the present problem, these fluctuating forces are represented by 
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the non-linear terms in the equation of motion due to the interaction of the system with the thermal bath when 
T ^ OK. The role of the fluctuating forces is to bring the system to the thermal equilibrium. The non-linearity of 
the equations, caused by the saturation of the atomic transitions, enhances the atom-bath interaction. 

Even though the works about fluorescence resonance take in account the interaction of the resonant system with a 
thermal bath, they in general assume the temperature of the bath as T — OK. Senitzky [T^] , however, in a elegant 
way regards the effect of thermal fluctuations in the resonance fluorescence of a two-level atom. But Senitzky treats a 
more general model for the thermal reservoir than we present here, and his mathematical treatment is more involved. 

In this paper we prefer to use the equations of motion of the dynamical operators to calculate the intensity of the 
fluorescent light. Although this approach is equivalent to the density operator, it has the advantage to facilitate the 
physical interpretation, because it resembles with the classical treatment given to the Brownian movement, in terms 
of the Langevin's equations. 

This paper is organized as follows. In Sec. 2 we describe formally the system and obtain the Heisenberg equations 
of motion in integral form. In Sec. 3 we determine the solution of the system of equations and apply it to the 
spontaneous emission to obtain the mean intensity of the scattered fields in the fluorescence of the two excited levels. 
In Sec. 4 we discuss and interpret graphically the results of the Sec. 3. Some details of the calculations are given in 
Appendix. 

II. THE HAMILTONIAN OF THE SYSTEM 

We consider here the problem of an atom fixed in space, with three levels and one electron, assuming that transitions 
occur only between each excited level and the fundamental one. Thus, the three-level system will be described with 
the help of the following operators 



Roo = |0)(0|, (2) 

Rn = |1)(1|, (3) 

R22 = |2)(2|, (4) 

R01 = |0)(1|, (5) 

R02 = |0)(2|, (6) 



which obey the relation 

1 = |0)(0| + |1)(1| + |2)(2|. 

The fundamental level is represented by |0), and the excited levels are represented by |1) and |2), respectively. The 
commutation relations of the Rij operators are 



[Roi,R 01 ] = Roo — R11, (7) 

[i?o2,-Ro2] = Roo — R22, (8) 

[i?02,i? 02 ] = -Rl2=0, (9) 

[RouRl 2 ] = -i?2i = 0, (10) 



where we make the transition rates between the two excited levels vanish. The atom is illuminated with two polarized 
laser beams; each beam has a frequency close to the characteristic frequency of each excited level of the atom. We also 
assume that the light beams are intense, and they will be treated classically. Within this point of view, the atom may 
be assumed as a couple of electric dipoles interacting with the electromagnets fields of the light beams. To account 
the radiative damping, a thermal bath is simulated with a big system of harmonic oscillators, that perform the role 
of the normal modes of the thermal radiation field. The Hamiltonian is 

H = Ho+H!, (11) 

where 



H = huoiRn + hw 2 R22 + fr/^^blbk. 

k 



(12) 
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b k and b\ are the operators corresponding to the modes of the bath (annihilation and creation, respectively), that 
satisfy 

[b k ,b{,] = 8 kk r. (13) 
The interaction Hamiltonian, in the rotating wave approximation, is given by 



H r = | — /ii4i UiEi^ + iJ^giityb^j -hRl 2 ^X 2 E' 2 (t) + i^2 g 2 (k) b k j^ + h.c. 



(14) 



where E[(t) and E 2 (t) are the driven fields, with frequencies respectively close to u\ and uj 2 . Xi and gi are coupling 
constants. The frequencies u>i and u> 2 are assumed quite different, so the driven field tuned with one frequency will 
not excite electrons to the level corresponding to the other one. The Heisenberg equations of motion of the operators 
are 



dR 



IK) 



tit 



dR n 
dt 

dR 22 
dt 

dR i 
dt 

dR 02 
dt 



-iR w ^Ai£;;(£) + z^.gi(fc)& fc j - 
-iR 20 ^\ 2 E' 2 (t) + i ^2 92(k)b k J | + h, 
iRw ^XiE[(t) + i^2 gi (k)b k ^ +h.t 
1R20 ^\2E 2 (t) + iY,92(k)bkJ I +h.( 
-iuj 1 R 01 + i{R oa - Ru) ^XiE[(t) +i^g 1 (k)b k 



(15) 
(16) 
(17) 

-mjx-Roi + i(Roo - Ru) I XiE[{t) + i} a^bk I (18) 



-ioj 2 



R02 + i(Roo - R22) \\ 2 E' 2 (t)+iY J 92ik)b k \ (19) 



^ = -iu' k b k +g* 1 (k)Roi+g* 2 (k)Ro2, (20) 

where h.c. means the hermitian conjugate. 
Let us define the Bi(t) operator as 

5iW = E5i( fc ) & ^), ( 21 ) 

and B 2 {t) as 

B 2 (t) = (k)b k (t). (22) 

k 

Integrating the equation ( |2(i| ) , we may write 

Bi(t) = B 01 (t)+ VI51WI 2 f : Roi(f dt + 

k J ° 

+ E 3i(%2* (fc) / Ro2(t')e^( t '- t )dt>, (23) 
k Ja 

B 2 {t) = B 02 (t)+ y^g 2 {k)gt{k) [ R 01 (t')e^( t '"') dt' + 
k J o 

+ El-9i( fc )l 2 f Ro2{t')e i ^i t '- t )dt', (24) 
Jo 
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where 

Boi(t) = ^gi(k)b h (q)e- i < t , (25) 

k 

and 

502^ = ^52^)^(0)6-^'. (26) 



We can now use the expressions ( |23D and fl24| ) in the equations (|l_5]) - (|19|), to eliminate the bk(t) variables of the 
thermal bath: 

= -iR 10 (XiE[{t) + iB^t)) - iR 20 (X 2 E' 2 (t) + iB 2 {t)) + h.c. (27) 

= {ii? 10 (A 1 E 1 (t)+iB 1 (t))} + / l .c. (28) 



dR 22 
dt 

dR 
~~dt 

dR 02 



= {iR 20 {\ 2 E' 2 {t)+iB 2 {t))} + h.c. (29) 



dt 



- = -iuxRot + i (i? 00 - Rn) (XiE[(t) + iB t {t)) , (30) 
-iuj 2 R Q2 + i {Roo ~ R22) {X 2 E' 2 {t) + iB 2 {t)) . (31) 



This is a set of non-linear differential equations that cannot be exactly solved. The non-linearity of these equations is 
due to the interaction between the three- level system and the radiation fields. They are the Langevin equations for 
the system and B\(t) and B 2 (t), by analogy with the Brownian movement, are the random fluctuating forces acting 
on the atom. To solve it we can try some approximation. As a first approximation, let us suppose that the interaction 
is sufficiently weak to be disregarded. Thus, the operators will evolve in time as 

Rm(t') = Rox(t)e- im ('-*), (32) 
R 2(t') = R 02 (t)e-^( t '- t ). (33) 

With these adiabatic approximations, the equations ( |23| ) and (|5|) becomes 

Bt(t) = B Q1 (t)+y2\ gi (k)\ 2 R Q1 (t) [ e i H-^)( t '- t )dt' + 
k Jo 

+ ^9i(k)9*2(k)Ro2(t) [ e^-^'-'W, (34) 
ft Jo 

B 2 (t) = B 02 (t) + y292(k)gi(k)R 01 (t) [ e ^-i) (*'-*) dt > + 
k Jo 

+ V \g 2 (k)\ 2 R 02 (t) I e^-^'-^dt'. (35) 
k J ° 



Assuming that the reservoir modes are very close, the sums in eqs. ( p4j ) and ( p5| ) may be substituted by integrations, 
where the number of modes in the frequency interval dui' is given by g(u)')du>' . With this consideration, and using the 
known result |H| 

1 V 

e l{u - u){t - t] dt' = -1— + tt6(lj' - u), (36) 



where V is the principal value of the integral, we may write the expressions B\(t) and B 2 (t) as 
Bx(t) = B 01 (t)+(ng(uj 1 )g 2 1 (u 1 )-iV J duo' 9 9 }^' ] ) R 01 (t) + 

\ J LO — L0 2 J 

B 2 {t) = B 02 (t) + (^g(uo 2 )g 2 2 (wa) - iV J ^'^^) i2oa(*) + 

+ Ufa)* maw - iv J dJ iM^M^\ Rm{t) _ 



Inserting now the expressions (gTj) and ( |38| ) in the equations of motion (16)-([L9|) we obtain 

dR u 



dt 



dR 2 ; 
dt 



dRo: 
~dt 



dR Q , 
dt 



iXiRwit)^ (t) - i\tR Q1 (t) E[* (t) - R w (t) B 01 (t) 
-flji (t)Ro! (t)-hR n (t), 

i\ 2 R 20 {t)E' 2 (t) - iX* 2 R Q2 (t) E' 2 * (t) - R 20 (t) B Q2 (t) 
-B\ 2 (t) R 02 (t) - k 2 R 22 (t) , 

-i(woi - iki/2)R 01 (t) - (fci2 (oj 2 ) - i5ui 12 ) R 02 (t) - 
+ {R n (t) - R 00 (i)) (Boi (i) - iX 1 E' 1 (<)) , 



-i(w 2 - ik 2 /2)R 02 (t) - (k 21 (ui) - iSuj 21 ) R 01 (t) 
+ {R22 (t) - Roo (*)) (-802 (t) - iA 2 ^ (t)) • 
In the above expressions we have used 



woi = — V / rfw 



or — u>i 



,p(u/)<? 2 2 (u/) 
W02 = uj 2 — V / aw 



ui — UJ 2 



y = ^P{^l)9x (^1), 



y = 7rp(w 2 )g| (w 2 ) , 
k 12 (uo 2 ) = up (u> 2 ) g x (uj 2 ) g* 2 (w 2 ) , 

&21 (^l) = I"/" (wi) .92 (^l) 5* (wi) , 



(5t^i2 (w 2 ) — V J duj 



,p{u o')gi (uj')g 2 (10') 
oj' — Ul 2 



' LU' — UJl 



The integral forms of the equations (|39|)-(|42|) are 

ft 

R n (t) = i?n(0) e - felt + i A 1 / i? 10 (t')^(^)e fel( *' _t) ^' 

Jo 

ft 

/ Rq! (*') E[* (t') e^'-^dt' - 
Jo 

- f R w (t')B Q1 (t')e k ^ t '~ t )dt' - 



o 



1 bI (t')R 01 (tV^-'W 



o 



JM*) = i? 22 (0) e- fc2t + zA 2 / R 20 (t')E' 2 (t')e k < t '- t )dt' 

Jo 

-*A* / i? 02 (*') Bj* (*') e^'-^dt' - 
Jo 

R20 (t')B 02 (t^e^'-^dt'- 







i?oi (<) = i? i(0)e- l ^ ol - lfcl / 2 ) t - 

- (^12 M - Ww ia ) / i?02 (*') e 4 ^ 1 - lfcl / 2 )( t '-*)dt / + 



(i?n («') - Roo (f)) (B01 (0 - *Ai^ («')) e^-*/^*'-*)^/ 



i? 02 (t) = fio2(0)e- ,( *- ,t!/2)t - 

- (fei (ui) - iScj 21 ) [ R m {t')e^- lk2/2) ^^)dt' + 



(R22 (f) - Roo (f)) (B02 (*') - *A 2 £; 2 (0) ei^-^/a)^-*)^ 



Supposing that the external fields are monochromatic and plane-polarized, we write 

£1 (t) = E m e-^\ 

E' 2 {t)=E m e- i "* t . 

Introducing this form of the external fields in eqs. ([5l|)-(|54]), and taking the mean values, we obtain 

(Ru (*)) = (Ru (0))e- felt + i A 1 £oi f (Rio ~*) dt - 

t 







iXlE^ / (i? i {t'))e- lUllt ' e^'-^dt' - 
Jo 

(R w (t')Boi (t'))e fcl ( t '- t )*'- 
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(R 2 2(t)) = (i?2 2 (0)}e- fc2t + i A 2 i?02 / {R2o(t'))e-^ t 'e k < t '- t )dt'- 

Jo 

-i\* 2 E* 2 / {R 02 {t'))e- lU2lt ' 'e^'-^dt' - 
Jo 

'{R20 (t')B Q2 (t')>e fe (*'- t )^'- 

(Bl 2 (t>)R 02 (t>))e k < t '- t )dt>, (58) 











<iZoi(<)> = (i?oi (0)} e -^ 01 - ,;fcl/2)t - 

-(k 12 {uj 2 )-i5u 12 {uo 2 )) [ (Ro 2 (t')}e i ^- ik ^( t '- t )dt' + 
Jo 

+ [ ((2R u (t') + R 22 (0 - 1) B 01 (*')) e^-oi-ifei/^Ct'-i)^/ + 
Jo 

-iAi^oi / ((2i?n (*') + i? 22 (*') - 1)) e-^'e^-^/W-Vdt', 



(59) 



(JZo 2 (*)> = (i?02 (0)}e- 4( ^- 4fc2/2)t - 

- (*ai (wi) - Wwai (wO) / (iJo!^'))^^ 02 "* 272 ^*'"*^*^ 
Jo 

+ / ((2i? 22 (t') + i? n (i') - 1) B 02 (t'))e i(wo2 - ife2/2) (*'- t )^ / + 
Jo 

-iX 2 E 02 / <(2J?22 (0 + R11 («') - l)} e - M2!t 'e j( " 02 ~ ife2/2) ( t '- t )dt'. 



(60) 



III. SPONTANEOUS EMISSION AND RESONANT FLUORESCENCE 



We follow the procedure outlined in ref. |l6| to obtain the solution of the system of integral equations ( |57}|6"c| ) . We 
use a iterative process: the expressions i?n, -R22, -R01 and i?02 are linear in the initial values of the operators and 
the coefficients of the operators are expressed in terms of powers of Bqi (t) and Bq 2 (t) . Supposing that the radiation 
field is in a stationary state and has a Gaussian distribution, the higher-order correlation function can be expressed 
in terms of the second-order onep6[, 

n 

(Bt(h) ■ ■ ■ Bl(t n )B (O ■ ■ ■ B (t[)) = J2U (BktjWi)) ■ (61) 

\3\ i=1 

The sum must be taken over all permutations j = 1, ... , n. As the expressions of the operators depend on the 
radiations field configuration in earlier times, the only terms that contribute to the series expansion of equations ( p37| ) 
and ( p^ ) are those for which the field operators are time ordered, 

min(tj,4)> max (t i+1 ,t' i+1 ) , i = 1,2, 

The only term in the sum of the equation ( |6l| ) that contributes is 

n 

1=1 

Thus, terms that contain a different number of factors B^(t) and Bo(t) do not contribute to the mean value of the 
operator. A further simplification can be obtained if we suppose that the radiation field spectrum is dense, flat and 
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broad. In this case, the correlation functions (_B{ (t")i?io(iQ) and (i?| (0-^2o(ii)) can be calculated, by using 



equations (|25[) and (p26|) , 

(BkOiMO) - E Isi Wl' * M )(*"-*'), (62) 

(5| (i")5 20 «.)) - E 1-92 Wl 2 N (w) e l( — -)('"-*'), (63) 

where 

7VM = (6t( )6 fc (0)> (64) 

is the initial mean number of photons at the k mode of the radiation field. Substituting the sums in equations (|6 
and ( [63] ) by integrations, we have 

(Bt Q (f)Bi Vi)) = J dupMsAWNMe^-^W'-n, (65) 

(^ (t")S 20 (^)) = rdup(w)gl(u>)N( U ) j<"*-**>W'- t '), (66) 



and using the properties of the Dirac delta- function and definitions (|45| ) and fl46|) , we obtain 

(B\ Q (t")B w (t')) = ^N^W - t'), (67) 



(Bl(t")B 20 (t')) = ^N(u> 2 )5(t" - t'). (68) 



Inserting eqs. (69) and the complex conjugate (|59| ) in eq. (p7|), we obtain 

(R n (t)) = (i?ii(o)) e - fci * + r 21 ( e -^*-e- feit ) + r; 1 ( e - 2 ^- e - feit ) + 

r 11 (l-e-' £lt )-2A :i 7V(w 1 ) / (i? 11 (i'))e' £l(t '- t) A'- 



-*iJV(wi) / (i?2 2 (t')}e' £l(t '- t) ^'- 

f^dt'e'^e^ 1 ' jjf' (i? 11 (t"))e Zlt "^" 

>2 ft ( ft 1 

(i? 11 (t"))e z i* di" 

^2 ft ( r t' 



~T J* dt ' e ~ klteZ ' lt ' \J* (W)) eZlt " dt " 



h 

2 

Now we can calculate the correlations that appear in eqs. ( f)7\ ) and (|58|). Using relations j6l|), (|67| ) and (jS8|), in 
expressions ( |5l| ) and (|52]), we obtain the expressions of (Rio(t')Boi(t')) and (R2Q{t')B m {t'))- 

{R w {t')B 0l (t')) = hiN{u3x){Rxi{H)) + y «i?22(i')> - 1) - (69) 

{R2 {t')B m {t')) = k 2 N{u 2 ){R 22 {t')) + ^ ({Rn(t')) - 1) . (70) 

Expressions (^) and ([70]) are real, and then, 

(R 10 (t')Boi(t')} = (BUt>)B m (t>)), (71) 

<i2 a o(*')Bo2(< , )> = <4» (t')Ro2 (*')>■ (72) 



* dlfe- klt e Zlt ' \ f (R22(t"))e z '^" dt"\ . (73) 



Inserting the equation ( |70| ) and the complex conjugate of eq. (60) in eq. (|58|), we have 
(R 2 2(t)) = (i? 22 (0)}e- fe2t + r 22 (e- Z2t - e- fc2t ) + T* 2 (e-^ - e~ 
T 12 (l^e-^ t )-2k 2 N(u J2 ) f (i? 22 (i'))e fc2(t '-* ) ^' - 



Jo 

~k 2 N{uo 2 ) [ (R^t'^e^'^df - 
Jo 

-^J^dt'e-^e^'lj 1 (R 22 {t"))e^ t " dt" 



O 2 /"* 



(i? 22 (i"))e Z2 * 4 di" 



o 2 /■* 

2 ' dt / e- fcat e^ v 



4 Jo 
O 2 /"* 



(i?n(i"))e Z2t di" 



<fln(t")>e* 2 *' dt" 



In expressions ([73) and (|74[) we use the below definitions 



r 11 = jv( Wl ) + 

4N 



r 12 = iv(^ 2 ) + 

4 2 2 



r 2 i - 



iAi£ i<-Rio(0)> n? 



4ki| 



r 22 



«A 2 ^ 02 (i? 20 (0)) ft 2 



4|z 2 



2 ' 



Hi = 2|Ai||^,i|, 

n 2 = 2|a 2 ||£ 02 |, 



z 2 = y + «(^ 2 ; - Woa). 



When the external fields vanishes, the equations ([73|) and ( |74| ) simplifies 

(i?nW) = ( J Rn(0)) e - fclt +iV( Wl )(f- e - fcl *) - 



kxN{u!) / ( J R 22 (i'))e fcl (*'- t )*', 
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(iM*)) - (R22(Q))e- k2t + N(u; 2 )(l-e' k - t ) - 

~2k 2 N{u; 2 ) [ (R22{t'))eM t '- t )dt / - 
Jo 

-k 2 N{uo 2 ) ( (i? 11 (t'))e fc2 ( i '- t )^'. (84) 



Integrating these equations, we obtain 



(R n (t)) = ( J Rii(0))e- felt + ^7V( a;i )(i_ e -'= 1 *) 



ki 

f ^N(uj 1 ){R 22 (0)) (e-^ - e^ 4 ) , (85) 



fci — k 2 

(R 22 (t)) = ( J R 22 (0)) e ^ 2t + ^N(u> 2 ) (l - e- l2t ) - 

-=-^J\TM<flii(0)> M 1 ' - e- 124 ) , (86) 



where 



*i = Ai (1 + 2JV(wi)) , (87) 
fc 2 = k 2 (1 + 2JV(wa)) • (88) 

If the initial state of the system is given by 

ip= |0)a + |l)ai + |2)a 2 (89) 

with |a | 2 + |ai| 2 + |a 2 | 2 = 1, then 

ki 



(R u (t)) = \ ai \ 2 e- k ^ + ^-N{Lo 1 )(l-e- k A 

k\ v / 



ki - h 



-N^^fe-^-e-^), (90) 



(R 2 2(t)) = \a 2 \ 2 e- k ^ + =^N(LU 2 )(l~e~ k A~ 

k 2 \ > 

-N(cj 2 )\ ai \ 2 (e^ -e~ k A. (91) 



k 2 - k\ 

If the temperature of the bath is zero, and the system is initially in one of the excited states, say, the state |1), then 
fci = kx, ciq = a% = and 

(R n (t)) = |ai| 2 e- fcl * (92) 

which is the expected spontaneous decay of the excited state with lifetime l/k\. If the temperature is different of 
zero, the thermal fluctuations enhance the coupling between the system an the field and the decay rate is increased, 
as we note by the equation (^0|). To times t 1/fci, the system approaches the saturation regime and then 



(93) 



Assuming that the radiation field is in thermal equilibrium with the cavity, we have 

JV(wi) ' 



e ft,wi/fe B T _ l 



and 



1 



and the equations (|90|) and fell) becomes 



(iZii(t = oo)) = 
(R 2 2(t = oo)) = 



1 



e hui/k B T _|_ ^ 
1 

e hu] 2 /k B T _|_ ^ 



which is the usual Fermi-Dirac distribution. 

With the procedure discussed in detail in Appendix, we obtain the solution of integral equations (173) and (p 



(i* u (t)> = (iJn(oo)> + (fl u (0)> 53 



Ufi + h/2) 2 + Au 2 



Jit 



+*iJV(wi) J! 



(/< + fci/2) 2 + Aw\ 



2 ,j,k— 



/i (/* - /i) (/< - A) 



zn 1 ( J R 10 (0))e^ 



(/« + 4) Jit 



(fi - fj) (fi - fk) 



(fi + z l) hi 



(/t _ /j) (/t _ fk) 



E 

i,j,fc=Q 

ir> 1 (fi 01 (0))e-^ A 
2 ^ 

ijtj^k 

JJ? y> (/i + fci/2) /4t _ 

3 + fcl / 2 ) 2 + Aw? (/i + fc 2 /2) 2 +Aw| 

-fciiV^i)^^)) ^ 



(/« - /j) (/i - fk)h(fi) 



-k 1 N(Lj 1 )(R 2 2(0)) 



3 + fc x /2) 2 + Aw? (hi + k 2 /2) 2 + Aujj 



Jit 



2 



f(hi) (hi ~ hj) (hi - hk) 



-O 2 <i? 22 (0)> 53 



3 (/. + fc 2 / 2 ) 2 + Awf (fi + ki/2) 



z,j,k=Q 



(fi-fj)(fi-fk)h(fi) 



-« 2 (i? 22 (o)) 53 



+ fc 2 /2) 2 + Aw 2 


(hi + kx/2) 




fi - fk) 



e fit_ 



,hit 
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and 



3 \(h l + k 2 /2) 2 + Aw| 

(R22(t)) = (i?22M) + (i? 22 (0)} ^ 



3 (ft, + fca/2) 3 + Au>§ 
+fc 2 iV(a; 2 ) ^ ^_ ^ _ ^ ^ _ ^ 

, iO 2 (i? 20 (0))e^ 



E 



i^j^k 



9! 



+ z 2 ) 



i,j,k=0 



2 3 



"2 



(fci + to/2) eM 



2 ^ ( h i ~ hj) {hi - h k ) 



3 \(hi + k 2 /2) 2 + Aw 2 (^ + fci/2) 2 + Aw 2 

- W( K ll( 0» _ £ J to-w-wM 

3 !(/. + fc 2 /2) 2 + Aw 2 (/, + A; 1 /2) 2 + A W 2 

**** <0,) , £ 1 M/.) (/.-/,)(/. -A) 



-^<i?n(0)> £ 



3 + fcx/2) 2 + Aw 2 (hi + k 2 /2) 



i,i,fc=0 



{hi - hj) {hi - hk) f{hi) 



-fil(i?u(0)> ^ 



(/ 4 + fc!/2) 2 4 


-Aw 2 


(/< + to/2) 









In the above equations, we have used Awi = u>u — u>\, /S.uj 2 — lu 2 i — lu 2 and 

(fc 2 /4 + Aw 2 )7V( Wl ) + f>2/ 4 



(-Rn(oo)) 



(fc 2 /4 + Ac; 2 ) (1 + 2iV(wi)) + flj/2 ' 



(99) 



(100) 



(fc 2 /4 + Aw 2 ) 7V( W2 ) + ^ 2 /4 
(fc 2 /4 + Aw 2 ) (1 + 27V(w 2 )) + 2 /2 ' 



(^M) = 7T2 //xA,.,2uLo^lo2 75 > ( 101 ) 



IV. DISCUSSION OF THE RESULTS 

In the equations (|98|) and (|99| ) the only terms that remains when i — > oo are (i?n(oo)) and (R 22 {oo)). The Fig. |l| 
shows the behavior of (i?n(oo)) as a function of the detuning Awi and of the occupation number N (w\). The 
system rapidly tends to the maximum intensity of scattered light for occupation numbers N (u>i) ^ 0, and the peak 
of resonance with the driving field becomes less effective to increasing number of occupation. 

We see by an exam of the equations ( p8| ) and (|9^) that one excited level is affected by the other only in transient 
terms that depend of the initial values of (Rn{t — 0)), (R 22 {t — 0)). Thus, within the second order approximation 
of the coupling constants we have used, the upper levels are almost independents, and the three-level system behaves 
as two two-level systems. Some numerical values given by equation ( |98[ ) for certain combinations of the parameters A 



13 



and 6 for the initial state (i?u(0)) = (i? 22 (0)} = and (-Roi(O)) = (-Rio(O)) = 0, are given by Figs. |- |. With these 
initial values the equation (|98| ) becomes 



3 (/j + ki/2) 2 + Auf 
(JZ u (t)) = (E u (oo)} + fc 1 AT(a; 1 ) ]T 



+ V (/ ' + fcl/2) e/«* (102) 

The Fig. |^ shows the effect of the A1.2 parameters, that are proportional to the intensities of the driving fields, 
over the fluorescent light intensities (Ru(t)) and (i? 2 2(^)) for Bose occupation numbers N(wi) = N(w2) = and 
9i=e 2 = 1. 

In the Fig. ||| the light intensities (Rn(t)) and (-R 22 (t)} are showed to the same parameters of the Fig. ||, but now 
with Q\ = 10 and 6> 2 = 3. It can be seen that the evolution in time of the fluorescent light intensity for N (wi) = is 
increasing oscillatory as the detuning increases. This feature is attenuated for N (u}\) ^ 0, due to saturation caused 
by the combined effect of the thermal bath and the driving fields. 

The Fig. || and Fig. || shows that a increase in the Bose occupation numbers N(w\) and -/V(w 2 ) increases too the 
light intensities, but the amplitude of the oscillations showed in Fig. || becomes relatively smaller. The system rapidly 
attains the thermal equilibrium with the bath, damping the oscillation amplitude. 

The intensity of the fluorescence light which for weak driven fields is determined by the temperature of the bath 
(Fig. |^) , for strong fields depends mainly of the intensity of the driven field (Fig. [?]) . 

By the preceding examples, we see that the presence of the thermal bath modifies strongly the time evolution of 
the fluorescent light intensity. Thus the spectral density and the intensity correlation of the fluorescent light must too 
be affected. We will present these calculations in a further paper. 

V. APPENDIX 

Equations ([73]) and ( |74| ) form a set of integral equations that can be solved by Laplace transformation, 

<«*(•)> = W) + r 21 W) + 

r* 1 (s + zi)(fc 1 -z*) fcirn(s + zi)(s + z?) 

/(*) */(*) 



[fciiV(wi) (s + zi) (s + z*) + Of (s + fci/2)] 



(103) 



,„ ( (fl22(0))(a + z 2 )(a + ^) , (ha ~ Z2) (s + z*) 

(i?22(s)) = hfa) +F22 W) 

(a + z 2 ) (fci - 4) fc 2 ri 2 (a + z 2 ) (s + z* 2 ) _ 
h(s) sh(s) 



h(s) 



z 2 ) (s + z|) + ^2 (s + fc 2 /2)] , 



(104) 



In equations (103) and (104) we have used 

f(s) = (s + fci) (s + Zi) (s + z*) + 2fci7V( Wl ) (s + zi) (s + z*) + n? (s + fci/2) , (105) 



h(s) = (s + k 2 ) (s + z 2 ) (s + z* 2 ) + 2k 2 N(uj 2 ) (s + z 2 ) (s + z|) + (s + fc 2 /2) 



(106) 
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Substituting (104) in (103), and retaining the terms until second order in the coupling constants, we have 



<i?ll(s)) = ?w + — + 

r^zi (s + zi) k%Tn (s + zi) (s + z^) _ 

/00 */(*) 

- (J?22(0)) / ( ( S |4;i (s+4) (s + *> (s + + ^ (s + fci/2) ] 

Taking now the inverse Laplace transform of the above equation, 

/■ioo+a 



-j pioo-\-a 

(Rn(t)) = — (R n (s))e st ds, t>0 

Z7TI J-ioo+Q, 



(107) 



where a is chosen in such a way that all singularities of the integrand fall in the left of the line Re s — a in the 
complex plane. Using the residue theorem, we obtain the wanted integral over the line Re s = a. Assuming that we 
may write 

/(*)= («-/i)0»-/2)0»-/3), (108) 

h(s) = (s-h 1 ){s-h 2 )(s-h 3 ), (109) 
where fa, fa, fa are the three roots of the cubic equation 

s 3 + [1 + N{ Wl )} 2k lS 2 + | [5 + 87V (wi)} +e 2 1 +X 1 \ f3fs+ 

+ |2 [1 + 2N ( Wl )} + 2 [1 + 2N («;!)] ft 2 + Ai| /3? = 0, (110) 
and foi, ft,2j ^3 are the three roots of the cubic equation 

s 3 + [1 + iV(w 2 )] 2k 2 s 2 + \[5 + 8N (w 2 )} + e 2 2 + X 2 \ f3 2 s+ 



(111) 
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where 0j = 2Auj l /k i , Aj = Of //3f , /3, = fcj/2 , with i = 1, 2. Then we have 



3 \(f i + fci/2) 2 + Aw 2 
<i2n(t)> - <flii(oo)) + <i2n(0)) (/i -/.)(/i-/ fc ) 



+kN(co 1 ) 



3 (/. + fcl / 2 ) 2 + Aco 2 



tni(flio(0))e^ A 
2 ^ 

»,j,fc=0 



(/i + *l) 



ifii(iioi(0)>e-« 



(/i - /;) (/i - A) 



(/» + zt) 



(fi - fi) (fi - fk) 



E 

i,j,k—0 

(fi + fci/2) /<t 



fi 2 

i=£j¥=k 

3 |( /t + fcl /2) 2 + A w 2 (/, + fc 2 /2) 2 + Aco 2 



-fclJV(wi)<iJ2 2 (0)> 



(/i, + fci/2) 2 + Aw 2 (hi + k 2 /2) 2 + Aw. 



i,j,k—0 
i^j^k 



f(hi) (h t - /ij) (/ij - ft fe ) 



-nl(R 22 (o)) 



3 (/s + fc 2 /2) 2 + Aw 2 (fi + fci/2) 



i,j,k— 



(fi- fi)(h- fk)h(fi) 



-n\(R 22 (Q)) J2 



3^ + fc 2 /2) 2 + Aw 2 (fti + fei/2) 



/(fc) (/i - £) (/* - /k) 
With a similar procedure we obtain a solution to (R 22 (t)). 



i,j,k=0 
i^j^k 



(112) 
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FIG. 1: (.Rii(oo)) dependence with the occupation number iV(aii), the detuning Acji and 



parameters f2i = 0.5 and fci = 1. 




FIG. 2: Time development of the fluorescent light intensities (Ru(t)) (points) and {R,22(t)) (solid line) for Bose occupation 
numbers N(wi) = N(w 2 ) = and Ai = 0.4, A 2 = 0.3, 9i = 62 = 1. 




FIG. 3: Time development of the fluorescent light intensities (Rn(t)} (points) and {R.22(t)) (solid line) for Bose occupation 
numbers N(wi) = N(w 2 ) = and Ai = 0.4, A 2 = 0.3, 6»i = 10 and 9 2 = 3. 




FIG. 4: Time development of the fluorescent light intensities {Ru(t)) (points) and {R,22(t)) (solid line) for Bose occupation 
numbers N(wi) = N(w 2 ) = 0.01 and Ai = 0.4, A 2 = 0.3, 9i = 10 and 6 2 = 3. 




FIG. 5: Time development of the fluorescent light intensities {Ru(t)) (points) and {R,22(t)) (solid line) for Bose occupation 
numbers N(wi) = N(w 2 ) = 0.1 and Ai = 0.4, A 2 = 0.3, 9i = 10 and 9 2 = 3 




FIG. 6: Time development of the fluorescent light intensity (Rn(t)) for Ai = 0.4, 9i — 10, and increasing occupation numbers 
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FIG. 7: Time development of the fluorescent light intensity (Rn(t)} for Ai = 10, 6 1 = 10, and increasing occupation numbers 



